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Abstract. - The dynamics of planar crack fronts in heterogeneous media is studied using 
a recently proposed stochastic equation of motion that takes into account nonlinear effects. 
The analysis is carried for a moving front in the quasi-static regime using the Self Consistent 
Expansion. A continuous dynamical phase transition between a flat phase and a dynamically 
rough phase, with a roughness exponent £ = 1/2, is found. The rough phase becomes possible 
due to the destabilization of the linear modes by the nonlinear terms. Taking into account 
the irreversibility of the crack propagation, we infer that the roughness exponent found in 
experiments might become history-dependent, and so our result gives a lower bound for £. 



Introduction. - The dynamics of cracks in heterogeneous media is a rich field encom- 
passing a large range of physical phenomena. In such situations a commonly studied quantity 
is the so called roughness exponent f. However, it is important to distinguish at least three 
different roughness exponents [1]: one describing the roughness in the direction perpendicular 
to the crack propagation, a second the roughness in the direction of the propagation, and a 
third one (which will interest us in the following) describing the in-plane roughness of the 
crack front during its propagation through the material. The exponent characterizing this 
in-plane roughness, has been measured in different materials, where it was found to be in 
the range 0.5-0.6, over at least two decades [2,3]. Despite numerous efforts [4-6], there is 
unfortunately no satisfactory theory that predicts the value of this exponent. 

In this paper we intend to contribute to the theoretical understanding of this problem. 
For that purpose we use an equation of motion of the crack front h(x) derived previously [7]. 
This equation contains two important ingredients - irreversibility of the propagation of the 
crack front and nonlinear effects. It is given by 

Qh 

— (x, t) = + h 12 [Kj (h) - K c (x, h)\ 6 [Kj - K c ] , (1) 



where h! — dh/dx, O(-) is the Heaviside function, Kj(h) is the stress intensity factor of the 
crack front (calculated to second order in h [7]) and K c {x, h) is a random term representing 
the heterogeneity in the local material toughness due to disorder. The random term can 
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always be separated as K c (x, h) = K* + rj (x. h), where K* is an average toughness and r) is 
its fluctuating part. It should be emphasized that terms including the effects of the system 
size were consistently left out in the derivation of this equation [7] . 

Solutions of stochastic growth models such as Eq. Q exhibit scaling behavior which is 
described using the time dependent height-height correlation function 



where C is the roughness exponent of the interface and z is the dynamic exponent. The 
brackets (• • •) denote average over disorder. In the following, we discuss several observations 
and arguments in order to explain our approach. This will help to simplify the equation of 
motion and to apply the self-consistent-expansion (SCE) approach in order to derive results 
for the scaling exponents. This method was developed by Schwartz and Edwards [8, 9] and 
has been applied successfully to the Kardar Parisi-Zhang (KPZ) equation [10]. The method 
gained much credit by being able to give sensible predictions for the KPZ scaling exponents 
in the strong-coupling phase above one dimension where many renormalization group (RG) 
approaches failed [11]. Another point which is especially relevant for our purpose is that for a 
family of models with long-range interactions (of the kind treated presently) SCE reproduced 
exact one-dimensional results while RG failed to do so [12]. 

The general approach. Regarding Eq. l|T)l. we expect to see three different regimes: 
A static regime for which K <C K* (where the Heaviside function in ||TJ can be safely 
approximated by 0); A regularly moving interface for large values of Kq (where the Heaviside 
function can be safely approximated by 1); And an intermediate complex regime, where 
Kq ~ K* . In this last regime, a very important factor seems to be the stabilizing terms which 
were dropped out in the derivation [7] , but which will make sure that the crack will stop after 
a while (as indeed seen in experiments [2,3]). 

Based on that picture, we hypothesize that a frozen dynamically rough interface is seen in 
experiments [3]), rather than a rough phase determined by a static pinned interface. In other 
words, we stress the point that the crack tends to stop due to its physical nature even without 
the presence of heterogeneities. This is indeed the case in cantilever beam experiments [3], 
where the crack faces are increasingly opened in order to induce crack front motion. As a 
result the front starts moving until it stops. The heterogeneities only induce roughness and as 
we argue, a dynamical roughness, which is then frozen due to the irreversibility of the fracture 
process. 

In order to test this picture, we approximate this system by neglecting consistently all 
mechanisms which deal with the slowing down of the interface, as well as the freezing of 
it. The assumption here is that the specific aspect of fine-tuning the opening stress mode 
(for example by imposing a time-dependent external loading) is exactly what the experiment 
does. Then we analyze the system at that critical point whichever means were taken to get 
there. This involves neglecting the Heaviside function on the right-hand side of Eq. JJ). We 
suspect that this term does play a role in the final stages of freezing, namely by imposing 
differential arrest along the interface (note again that the interface would stop anyway, even 
without this term). This would tend to increase the roughness. Thus, we would consider the 
results obtained below as a lower bound for the roughness, offering a quantitative physical 
explanation up to the last steps of the freezing. 

Interestingly, a similar approximation is implicitly present in the KPZ system. It is well 
known that any rough surface would eventually flatten by the KPZ system if the noise is 
stopped [13]. However in real situations, it is compensated by non-zero "angle of repose" 
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that eventually freezes the system in a rough phase [14] (this is expressed by an additional 
Heaviside function in the KPZ equation). It is also shown that the roughness exponent would 
be the same as that of the driven system if the freezing is done adiabatically [14]. This shows 
that a hysteretic effect (existence of an angle of repose) can be consistently neglected once an 
above-threshold driving noise is present. It also hints that the final roughness of our system 
would be "history dependent" i.e. it would depend on the protocol of the loading/freezing, if 
not done adiabatically. 

Following the previous arguments, we approximate the noise term for the moving front, 
where h ~ vt, by T)(x,h) ~ rj{x,vt) — fj(x,t) [15]. Also, we do keep nonlinear terms, since 
we claim (and will justify later) that they play an important role in roughening the interface. 
Obviously a linear equation of the kind described above (i.e. taking into consideration only 
the linear term in Kj(h)) would not yield any roughness, and actually even if the KPZ 
nonlinearities (i.e. h terms) are kept, we would also end up with a smooth surface, or at 
most logarithmically rough (this is a special case of the so called Fractal KPZ equation studied 
previously in [16]). When keeping consistently second order terms, the resulting equation of 
motion becomes 



dh f h'(x') dx' f dx' f dx" ti (x')ti (x") 

— = K Q , , v \ — +K ' 1 



dt " J (x' - x) 2tt u J 2vr J 2tt (a;' - x) (x" - x') 

— o> 

3 fA 



— oo -co 



s K*-K \h' 2 + (K -K*)+fj(x,t) , (3) 

with noise correlations described by 

(?) (z, t) f) [z\ t')) = 2D5 (z - z 1 ) S(t- t') , (4) 

where D is the variance of the noise. The integrals in Eq. J2| should be taken in the sense of 
Cauchy principal value. The first prediction of this equation is that the mean velocity should 
be proportional to the constant term, i.e. v cx (Kq — K*). The natural appearance of a velocity 
strengthens the simplification of taking a time-dependent noise term. For convenience, the 
constant term can be scaled out by transforming into a co-moving coordinate system, i.e. 
h — * h + (Kq — K*)t. Then, by looking at the KPZ term (i.e. h' 2 ) we can estimate the 
region where this discussion is relevant. Roughly, when the coefficient of that term remains 
negative (i.e. for Kq < |if*), we are still in the quasi-static regime since in that case a rough 
interface would decrease the velocity, while for higher values of the applied stress (-Ko > ^K*) 
the system would be in the regularly moving regime. This estimate is consistent with our 
assumption that the dynamics of interest is not necessarily at Kq ~ K* , but in some range 
above it (i.e. K* < Ko < ^K*). In the following, we will neglect the KPZ-term to simplify the 
presentation. We checked that the results we get for the scaling exponents do not depend on 
its existence (beyond its stabilizing effect). At the end, we will comment on the modifications 
due to its presence. 

The SCE method. - The SCE method is based on going over from the Fourier transform 
of the equation in Langevin form to a Fokker-Planck form and on constructing a self-consistent 
expansion of the distribution of the field concerned. We then consider the simplified version 
of the equation of motion in Fourier components 

dh q (t) _ , 



dt ~ q q ^ 



hg-y. M q i m hih m + fj q (t) , (5) 
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where c„ = \q\ and M q t m — — \q\ \i\ Sq.e+ m , L being the size of the front. Note that 



in contrast to the KPZ problem M q i m has the symmetries M q g m = M^ q j^ m = M q _i 
Mg^-m. Last, fj q {t) is a noise term with zero average described by its variance (f) q (t) f) q i (£')) 
2DS (q + q') S (t — t'). Rewriting this equation in a Fokker-Planck form we get 



dP_ >p d 
dt + ^dh q 



A)T^ 1" C q h q + ^ M qtm h t h r , 



P = , (6) 



where P ({h q } , t) is the probability functional for having a height configuration {h q } at time 
t. 

The expansion is formulated in terms of the steady-state structure factor (f> q = (h^ q h q ) (or 
two-point function), and its corresponding steady-state decay rate that describes the rate of de- 
cay of a disturbance of wave vector q in steady state, namely w^ 1 = J °° (h- q (0)h q (t)) dt/ (h- q h q ). 
From the scaling form (0) it follows that for small q's <p q and uj q behave as power laws in g, 
namely <\> q = A|g|~ r and uj q = B\q\ z , where z is the dynamic exponent, and the exponent T 
is related to the roughness exponent by ( — (T — l)/2. 

The main idea of SCE is to write the Fokker-Planck equation dP/dt = OP in the 
form dP/dt = [O + 0\ + O2] P, where O , 0\ and O2 are zero, first and second order 
operators in some parameter. The evolution operator Oq is chosen to have a simple form 

Oq = — gf- h ] , where D q /oj q — cf> q . Note that <p q and u> q are still un- 

q « V -9 / 

known. Next, an equation for the two-point function is obtained. The expansion has the form 

4> q — 4> q + e q {<f)p,uip}, where e q is a functional of all </>'s and w's. This reflects the fact that 

the lowest order in the expansion is exactly the unknown <f> q . In the same way, an expansion 

for u) q is given by u> q — uj q + d q {(f> p ,uj p }. Now, the two-point function and the characteristic 

frequency are determined by setting e q {4> p ,ujp} — and d q {4> p ,ujp} = 0. To second order in 

the expansion, we get the following two coupled integral equations 

A, - ^-\q\ct>q+h{q)<t>q + h{q)=Q, (7) 



^ - ^\9\ + J(q)=0, (8) 



with 



h (q) = H M J dl |£| KI(k "' + '^^:: l(m+M) ^ . ( 9 ) 

j(o) = K ° \ a \ f d£\f\ W\i- i \+\i\)'f> q -e+\q-^(W+k\)<l>e nil 

327T / u> e +u> q _ e ' V J 

It is interesting to mention here that Eq. Q can be understood as emanating from the short 
time balance of the original equation, while Eq. (jSJ comes from its long time balance [9]. 

These equations can be solved exactly in the asymptotic limit (i.e. for small g's) to yield 
the required scaling exponents governing the steady-state behavior and the time evolution. 
The difficulty here arises from the fact that the integrals involved, Ii(q), I%{q), and J(q), 
have contributions from large £'s as well as from small i's. Therefore, one must consider 
the contribution of the large I integration domain on the small q behavior of the integrals 
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<|9I11[) . For this, we break up the integrals Ii(q) and J(q) into the sum of two contributions 
if (q) + if (q) , and J < (q) + J > (q) , corresponding to domains of I integration with low and 
high momentum respectively. We expand If (q) and J > (g) for small q's and obtain the leading 
small-g behavior of the integrals 

I>(q) = |g| Igl^-dQ 2 , (12) 

If (<?) - A 2 q 2 - B 2 q 2 LU q + C 2 \q\ 3 , (13) 
J> (q) = A 3 \q\ - B 3 q 2 , (14) 

where the coefficients generally depend on the cutoff. Note that the constants A\, A 2 , A3, B\ 
and B 2 are strictly positive. Using these results, Eqs <|7I8[) reduce to 

9 _ ( KO _ . \ I , , , 7</„W , T< 



Do + A 2 q 2 A 1 \q\ <P q + If (q) <t> q + If (q) = , (15) 

!y-A^\q\ + J<(q) = 0. (16) 

The advantage of Eqs I|15ll6f) over Eqs (|7I8I) is that at the mere price of renormalizing some 
constants in both equations, we are left with the integrals If(q), I 2 (q) and J < {q) that can 
be calculated explicitly for small q's since the power-law form for (f>£ and tug for small ^'s can 
be used. The treatment of Eqs I|15ll6fl is carried on by studying the various possibilities of 
balancing the dominant order for small q. Note also that the small g-dependence of each of 
the integrals lf(q) and J < {q) depend on the convergence of the integrals without cutoff. So, 
to leading order in q we get 



^(9) ~ ( Ffqf-r-* 4-r-z<0 (17) 

T<(a) ~ / E2q2 3-2r-z>0 

2 {q) { F 2 \q\ 5 ~ 2r - Z 3-2r-z<0 (18) 

7<(a) ~ i EM 4-T-z>0 



We consider now the quadrant of the (r, z)-plane defined by Y > and z > 0, where 
solutions may be expected. The lines 4 — Y — z = and 3 — 2r — z = divide this quadrant 
into four sectors. The classical way [8] is to investigate each sector separately to decide whether 
or not a solution might exist there. After performing this analysis we can show that solutions 
are possible only in the sector defined by 4 — Y — z > and 3 — 2Y — z < 0, and so we present 
a detailed analysis for that sector only, where Eqs i|15lltjfl are rewritten as 

Do + A 2 q 2 -a(^--A 1 - Ek) M 1 ^ - BBi \q\ 1+z - r + C\ M 2 ~ r + F 2 q 5 - 2r - z = (Bp) 

B\q\ z - - A 3 - E 3 ^j \q\ - B 3 \q\ 2 = . (21) 

From Eq. I|21() , it can be easily seen that (in the small q limit) z = 1 corresponds to a possible 
solution. Then using this result in Eq. H2()|l. we find that either r = 1, or 1— Y = 5 — 2Y—Z. The 
last option implies Y — 3 corresponding to a roughness exponent £ = 1, which is inconsistent 
with the assumption of small gradients used to derive the equation of motion. Therefore, we 
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are left with T = z = 1, which is consistent with the defining condition of this sector, and is 
just the solution of the linearized equation of motion derived from (J3J). 

Interestingly, a more careful inspection reveals another option ignored at first sight, namely 
that of getting a fine-tuned case where the renormalized nonlocal elasticity vanishes. In 
Eq. ((201 or HU, this corresponds to the case K /2 - A x - E\ =0 or K /2 - A 3 - E 3 = 
respectively (note that the sign of the coefficients is crucial for that argument - and this is a 
particular feature of this system which does not happen in the KPZ system [8,10] for example). 
This fine-tuning corresponds to being on the critical point. 

Repeating the analysis with these options in mind gives rise to three new possible phases. 
First, it is possible that the coefficient of |(?| 1_r vanishes in Eq. H20JI . implying T — 2. Then 
Eq. (|21[) implies z = 1 or z = 2 according to the coefficient of |g|. Another option is when 
only the coefficient of |g| vanishes in Eq. (|21() . which implies z = 2, while in Eq. (|20|l as before 
r = 1. A more careful estimate shows that A x + A 3 < E x + E 3 and F 2 > 0, and so the 
only possible cases are with z = 2. In addition, since T = 2 yields a more singular balance 
in Eq. I|21|l (in the sense that the leading power of q becomes —1) compared to the one with 
r = 1, we expect that this scenario will be the dominant one. 

Summarizing this part, we found two possible phases: First, a flat phase described by 
C — and z = 1, corresponding to the system in the moving regime. This phase is always 
possible. Second, we see the possibility of having a rough phase with £ = 1/2 and z = 2, 
which is possible only on the critical point. 

Discussion. - In this paper we analyzed a recently proposed equation of motion for an 
in-plane crack front with the aim of studying possible roughening of the front. We found 
the possibility of having a rough moving phase with £ = 1/2 (and z = 2) which is relevant 
for Ko ~ K* due to destabilization of the nonlocal elasticity by the nonlinear term. This 
result is in agreement with the roughness exponent measured in experimental systems [2,3]. 
Since in our analysis we neglected the irreversibility of the fracture process (which becomes 
important during the last steps of freezing, and so tends to further roughen the line), our 
analysis provides a lower bound for the experimental results (recall that experimental results 
vary between 0.5-0.6). We hope that analysis of the full equation would yield results which 
are even closer to the experimental measurements. 

At this point it is useful to comment on what would have happened if we applied the self- 
consistent expansion to the full equation of motion © including the KPZ-terms. We checked 
and found that the basic analysis is not modified, with the only difference that the option 
of having a rough phase with £ = 1/2 and z = 1 is not ruled out like in the simplified case 
we discussed above (due to the fact that some prefactors change). This is interesting in view 
of that fact that z has been measured experimentally in [17] and found to be z = 1.2, close 
to the predicted value. This suggests that the KPZ nonlinearity plays an important role in 
the dynamics, and cannot be neglected. Finally, we think that these results are robust in the 
sense that they are applicable to other systems having a similar structure. 

* * * 
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